Incompressible strips in dissipative Hall bars as origin of quantized Hall plateaus 
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We study the current and charge distribution in a two dimensional electron system, under the 
conditions of the integer quantized Hall effect, on the basis of a quasi-local transport model, that 
includes non-linear screening effects on the conductivity via the self-consistently calculated density 
profile. The existence of "incompressible strips" with integer Landau level filling factor is investi- 
gated within a Hartree-type approximation, and non-local effects on the conductivity along those 
strips are simulated by a suitable averaging procedure. This allows us to calculate the Hall and 
the longitudinal resistance as continuous functions of the magnetic field B, with plateaus of finite 
widths and the well-known, exactly quantized values. We emphasize the close relation between 
these plateaus and the existence of incompressible strips, and we show that for B values within 
these plateaus the potential variation across the Hall bar is very different from that for B values 
between adjacent plateaus, in agreement with recent experiments. 

PACS numbers: 73.20.-r, 73.50.Jt, 71.70.Di 



I. INTRODUCTION 

The application of the quantized Hall effect (QHE)^ 
as resistance standard, and its importance for modern 
metrology, relies on the extremely high reproducibil- 
ity (better than 10~^) of certain quantized resistance 
values^ This extreme reproducibility points to an uni- 
versal origin, which is independent on special material or 
sample properties. The purpose of the present paper is to 
propose and evaluate a quasi-local transport model that 
allows us to calculate, first, the potential and current dis- 
tribution in a two-dimensional electron system (2D ES) 
under the conditions of the QHE, and, second, the longi- 
tudinal and the Hall resistance, Ri{B) and Rh{B), in the 
plateau regimes of the QHE and in between. Whereas 
the resistance values between the QH plateaus will de- 
pend on details of the used conductivity model, the ex- 
actly quantized plateau values result from the existence 
of sufficiently wide "incompressible strips" along which 
the local conductivity vanishes, since occupied and un- 
occupied states are separated by an energy gap (Landau 
quantization). Localization assumptions j^ which played 
an important role in early approaches to the QHE, are 
not included in our model. Our model is motivated by 
a recent experimental investigation of the Hall-potential 
in a narrow Hall bar^^ and a critical reexamination of a 
subsequent model calculation.'^ 

Experimental information about the actual current 
and potential distribution in a Hall bar under QHE 
conditions has been obtained recently by Ahlswede and 
coworker a'^'TiT i with a scanning force microscope^. The 
data were interpreted in terms of "incompressible strips" 
with constant electron density (corresponding to the fill- 
ing of an integer number of Landau levels) ^SiiSiiiii? which 
are expected to develop in an inhomogeneous 2D ES as 
a consequence of its strongly non-linear low-temperature 
screening propertiesii^ in a strong perpendicular magnetic 
field. If the filling factor in the center of the sample was 
slightly larger than an integer, the Hall potential was 



found to drop completely across two strips, while being 
constant elsewhere. With decreasing B, the strips moved 
towards the sample edges, just as one expects for the in- 
compressible strips in a sample, in which the electron 
density decreases gradually from a maximum in the cen- 
ter towards the edges. If the center filling factor was 
slightly below an integer, a gradual potential variation 
was observed, either linear or a non-linear, without clear 
indication for incompressible strips^ 

This interpretation was supported by subsequent the- 
oretical work of Giiven and Gerhardts (GG)5i, who 
extended the self-consistent Thomas-Fermi-Poisson ap- 
proximation (TFPA)ii*i2iii for the calculation of elec- 
tron density profile and electrostatic potential to a non- 
equilibrium situation with a position-dependent electro- 
chemical potential, determined by the presence of an ap- 
plied dissipative current through the sample. Electro- 
chemical potential and current density was calculated 
from a local version of Ohm's law, with a local model 
for the conductivity tensor, determined by the position- 
dependent electron density. The feed-back of the current 
distribution on electron density and the measurable po- 
tential profile^ was included by the assumption of local 
equilibrium in the stationary non-equilibrium situation. 

In agreement with the experiment **, the calculation^ 
shows a linear variation of the Hall potential across the 
sample if there are no incompressible strips, e.g. for suf- 
ficiently high temperature or if the magnetic field is so 
strong, that the local filling factor is everywhere in the 
sample less than two (spin-degeneracy is assumed and 
interactions which might lead to the fractional quantized 
Hall effect are neglected). Also for center filling factors 
slightly larger than 2 or 4 the calculation confirms the ex- 
periment, showing that the potential drops across broad 
incompressible strips and is constant elsewhere. How- 
ever, due to the use of the TFPA, GG- obtain incom- 
pressible strips whenever the center filling factor is larger 
than 2;2iiSiiiiiS and due to the strictly local conductivity 
model these dominate the current and potential distribu- 



tion, and lead to vanishing longitudinal resistance. Thus, 
the model assumptions of Ref 5 lead to serious disagree- 
ment with important aspects of the experiment. 

The purpose of the present paper is to improve on the 
model of GG^ so that, first, qualitative agreement be- 
tween the calculated and the measured potential distribu- 
tion is achieved for all filling factor regimes, and, second, 
reasonable results for Ri{B) and Rh{B) are obtained. 
Following the lines suggested by GG^, we investigate in 
Sec. m the conditions for the existence of incompress- 
ible strips, using a Hartree approximation. In Sec. IIIII 
we reexamine and weaken the strictly local conductivity 
model, and show that a simple spatial-averaging proce- 
dure of the local conductivities can simulate corrections 
expected from a Hartree calculation for the equilibrium 
state and from a non-local transport calculation. Trans- 
port results based on the self-consistent Born approxi- 
mation will be presented and discussed in Sec. IIVI In 
the present work we will restrict our consideration on 
the linear response regime, so that heating effects, which 
might destroy incompressible strips in the presence of 
high currents^ can be neglected. 



II. EXISTENCE OF INCOMPRESSIBLE 
STRIPS 

A. Electrostatic self-consistency 

Following Ref|a, we consider a 2D ES in the plane 
z — Q^ with translation invariance in the y direction 
and an electron density nc\{x) confined to the interval 
—d < X < d. The confinement potential Vhg{x) is deter- 
mined by fixed background charges and boundary con- 
ditions on metallic gates. The mutual Coulomb interac- 
tions between the electrons are treated in a Hartree-type 
approximation, i.e., are replaced by a potential energy 
term Vh{x) which is determined via Poisson's equation 
by the electron density. Exchange and correlation effects 
are neglected, and spin degeneracy is assumed. Thus, the 
electrons move in an effective potential 



V{x)=Vbg{x) + VH{x) 

K 



(1) 



Vnix) 



dx' K(x, x')nei{x'), 



where — e is the charge of an electron, R, an average back- 
ground dielectric constant, and the kernel K(x, x') solves 
Poisson's equation under the given boundary conditions. 
Kernel and background potential for the frequently used 
j^Q^g^^9ji0jiiji2 -j-j^g^^ assumes all charges and gates to be 
confined to the plane = 0, are taken from Ref \ 



where eno is the homogeneous density of positive back- 
ground charges in the Hall bar. Other meaningful and 
tractable boundary conditions are also possible. ^^ 

To perform explicit calculations, one needs a prescrip- 
tion to calculate the electron density for given effective 
potential V{x), which then together with Eqs. (Q) and 
^ completes the electrostatic self-consistency. The self- 
consistent TFPA^iiiiiSii^ takes this prescription from the 
Thomas-Fermi approximation (TEA) 



nei(a;) = j dE D{E) f{E + V{x) - fi*), 



(5) 



with D{E) the density of states (DOS), /(£') = 
l/[exp{E/kBT) + 1] the Fermi function, /i* the elec- 
trochemical potential (being constant in the equilibrium 
state), ks the Boltzmann constant, and T the tempera- 
ture. 



B. Hartree approximation 

A less restrictive approximation would be to insert 
V{x) into Schrodinger's equation, 



2m V c 



V{x) 



$A(r) = Ex^xir), (6) 



with A(r) a vector potential describing the magnetic field 
B = (0, 0, B) = V X A, and to calculate the density from 
the eigen- energies E\ and -functions $A(r), 



nci 



ir)=J2\^,ir)\'fiE,-^n- 



(7) 



Exploiting the symmetry of our system, we may use the 
Landau gauge, A(r) = (0, Bx, 0), and factorize the wave- 

functions, $A(r) — Ly eyip{iky)(j)n,x{x), with Ly a 
norma lization length, X = —Pk a center coordinate, 
/ — y^h/mojc the magnetic length, and ujc — eB/{mc) 
the cyclotron frequency. The Schrodingcr equation then 
reduces to 



'ff^ + '^c.lix-Xr + Vix) 
2m dx'^ 2 



4>n,x{x) = 

E,,{X)4>n,x{x), (8) 



(2) and the electron density becomes 



nc\{x) 



9s 
2ttI'' 



J2 J dXf{E„{X) - fi*Mr.,x{x)\\ (9) 



where gs — 2 takes the spin degeneracy into account 
and the sum over X has been replaced by an integral, 
Ly'Ex^i^nPr'JdX. 



K{x,x') 



In 



v/(d2-a;2)(d2_a;'2) _^ ^2 _ ^/^ 



{x — x')d 

^P - 



(3) 



Vig{x) = -ElJl-{x/df, El^27re^nod/R, (4) 



C. Thomas-Fermi Approximation (TFA) 

If the potential V{x) varies slowly on the scale of the 
magnetic length /, its effect on the lowest Landau levels 



(LLs) can be treated perturbatively, with the lowest order 
result 

E„{X)^E,, + ViX), E„^hLu,in + l/2). (10) 

On the length scale relevant for the variation of V{x), the 
extent of the Landau wavefunctions may be neglected, 
\4'n,x{x)\'^ ~ S{x — X). Then the Hartree result for the 
electron density, Eq. ^, reduces to the TFA, Eq. JSJ, 
with the Landau DOS 



D{E) 



Til 



^ oo 



En). 



(11) 



n=0 



To evaluate the self-consistent TFPA we follow Ref|3. 
First we fix the sample width 2d and the density of pos- 
itive background charges no, and thereby according to 
Eq. Q) the background potential Vbg (x) and the relevant 
screening parameter asc = T^ao/d, with oq = Rfi^ / (2me^) 
the screening length (2ao = 9.79 nm for GaAs). Next we 
choose the actual width 2b of the density profile at T = 
and B = 0, and solve for |a:| < b the linear integral 
equation^ 



V{x)-Vb,{x)^- 



dx' 



Kix,x')[^il^-Vix% (12) 



[with ijq = V{—b) — V{b)] to which the self-consistent 
TFPA reduces in this limit. From the corresponding 
density profile nci{x; B = 0,T — 0) = -Do[Mo ^ ^(2^)]; 
with Do = m/{TTh^) the DOS of the 2D ES at B = 0, 
we calculate V{x) for |x| < d, the average density 
fic\ = J_^ dx nci{x; B — 0, T — 0)/2d with the correspond- 
ing Fermi energy Ep = Uci/Dq^ and, for later reference, 
the Fermi energy E'^ = nci{0;B = 0,T = 0)/Dq corre- 
sponding to the electron density at the center. 

In the following we will consider only symmetric den- 
sity profiles and take &, or equivalently the depletion 
length d — 6, as a free parameter, that fixes the density 
profile and the electro-chemical potential fj,Q a.t B — 
(where the temperature dependence is weak). In real 
samples /ij ruay be determined by an electron exchange 
between the 2D ES and its surrounding, which may be 
possible at high but not at low temperatures. A restric- 
tion that fixes Hq will also determine the value of b. 

Next, we fix the value of the magnetic field and start 
with a high temperature to calculate the electron den- 
sity from Eqs. Q and (|ll|l self-consistently with Eq. ^, 
using the previously calculated potential V{x) as initial 
value. Finally we lower T stepwise until the required 
low temperature is reached, and iterate (using a Newton- 
Raphson approach) at each temperature until conver- 
gence is achieved. 

The solid lines of Fig. Q] show results for d — 1.5 /im 
and no = 4 • 10" cm"^ (which implies E^^ = 4.38 eV) 
obtained for 501 equidistant mesh points, —d — xq < 
x„ < xjq = d [N = 500). The density profile was fixed 
by choosing b/d = 0.9, which yields fiei = 2.9 • 10^^ cm~^ 
andrioi(0;B = 0,T = 0) = 3.58- 10" cm'^, and thus, with 



3 









2.02 














TFA 












/ / 






QHA 








/ 


■ / 






FHA 








/ / 


/ 










2 






J 


- 




^- — 




/ ,->'' 


y/' 




( 


/ 






/ 






/ 






I 


1.98 
-0 


\ 








85 


-0.84 -0.83 


Yo.82 


1 ^ 




"~-^ 


1 / 2.01 










\ 






' 






/ 






/ 




\ 




- 


/ ■' 






.^y 






- 


/' 








1.99 


/ 












.5 


-0.4 


-0.3 









-1 



-0.5 





x/d 



0.5 



FIG. 1; Electron density profiles for two values of the mag- 
netic field (hwajEp = 0.94 and 0.65) and different approxi- 
mations: Thomas- Fermi (solid lines), Hartree (dashed), and 
quasi-Hartree (dash-dotted). The insets show the enlarged 
plateau regions for both cases, a = 0.01, ksT / E% — 0.0124. 
Kinks in the upper inset indicate mesh points. 



Do = 2.8 • 10i°meV"icm-2 for AlGa, Ep = 10.37 meV 
and E^ — 12.75 meV. We prefer to use E^ (rather than 
E^ ) as an reference energy, since it has the same order 
of magnitude as the cyclotron energies of interest. For 
finite B we define an effective center filling factor vq = 
2TTpneiix = 0; B = 0,T = 0) = 2E°p/fuVc. 

The result obtained for the TFPA (solid lines in Fig.[TJ 
shows the expected well developed incompressible strips 
with constant electron density at local filling factor 
j^(x) = 2. For the larger B value we obtain wide den- 
sity plateaus at 0.32 < \x\/d < 0.46, in each of which 
we find at 36 Xn values, with high accuracy, ^{xn) = 2. 
For the lower B value the incompressible strips are much 
narrower, however we obtain the high precision values 
t^{xn) ~ 2 still at five neighboring Xn values. Typical re- 
sults of the TFPA are summarized in Fig.|21 which shows, 
as a gray scale plot, the filling factor profile for varying 
magnetic field, with horizontal lines corresponding to a 
fixed B value. At sufficiently large B field, the local fill- 
ing factor i^(a;) is everywhere in the Hall bar less than 2, 
and the 2D ES is completely compressible. At somewhat 
lower B (fkOc/Ep « 1) the center of the sample becomes 
incompressible with local filling factor i'{x) = 2, while 
I'ix) gradually decreases outside the incompressible cen- 
ter and falls off to zero in the depletion regions at the 
sample edges. With further decreasing B, the filling fac- 
tor in the center increases and incompressible strips with 
vix) = 2 move towards the sample edges and become 
narrower. At sufficiently low B values, incompressible 
strips with local filling factor i/{x) — 4 occur, first in the 
center, and then move towards the edges. They then co- 
exist with incompressible strips of i^(cc) = 2, which exist 
near the edges and are narrow, but, within the TFPA, 
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FIG. 2: (color online) Gray scale plot of filling factor ver- 
sus position X and magnetic field B, calculated within the 
TFPA. The regions of incompressible strips with i'{x) = 2 
and v{x) = 4 are indicated. For sufficiently large B {Qc = 
tkOc > E%) the system is compressible (indicated by "CS"), 
while for the lower B values included in the figure it always 
contains incompressible strips ("IS"). The dashed horizontal 
line refers to Fig.^lbelow; a = O.Of, kBT/E% = 0.0f24. 



still have a finite width. For low enough temperature, 
this type of behavior continues at stiU lower values of 
B, where further incompressible strips with successively 
higher fiUing factors evolve from the center and move 
towards the edges, coexisting with the edge-near incom- 
pressible strips of lower local filling factors. 



D. "Quasi- Hartree" Approximation 

The dashed lines in Fig.^are calculated in the Hartree 
approximation. We started again at B = 0, T = 
and inserted in each of the following iteration steps the 
previously calculated potential V{x) into the eigenvalue 
problem of Eq. (|SJ), took each mesh point a;„ as center 
coordinate X and diagonalized the problem in the space 
spanned by the eight lowest-energy unperturbed Landau- 
Hermite functions 



€A^) 



exp(- 



Xf/2f 



^JT^n\^l 



■B^{{x~X\H), (13) 



where H„(^) is the n-th order Hermite polynomial. The 
resulting energy eigenvalues and -functions were used to 
calculate the electron density according to Eq. 0. The 
overall appearance of the Hartree results for the density 
profiles and also the wide plateaus for the higher B value 
(see lower inset of Fig. ^ are in good agreement with 
those of the TFA. The narrow plateaus, obtained in the 
TFA for the lower B value, are now however smeared out. 
As is clearly seen in the upper inset of Fig.^ the Hartree 
result for the filling factor v(x) crosses the value 2 with 
a finite slope. 



FIG. 3: Magnetic-field dependence of the width of the v = 
2 incompressible strips, for three different approximations: 
the analytical result of Ref|i (CSG), the TFPA (TFA), and 
the quasi-Hartree approximation (QHA). Note the inverted 
B scale. Inset: ratio of the incompressible strip width to the 
magnetic length in QHA. 



The essential qualitative difference between the 
Hartree approximation and the TFA is the neglect of the 
extent of the wavefunctions in the latter. So we interpret 
the smearing-out of narrow incompressible strips in the 
Hartree approximation as being due to the finite width 
of the wave functions. To check this idea, we consider a 
"quasi-Hartree" approximation (QHA) in which, instead 
of solving the problem of Eq. (jSJ, we replace the wave- 
functions by the eigenfunctions 113|l of the unconfined 
Landau problem and take the energy eigenvalues from 
Eq. (|10|l . The latter would be correct in the sense of a 
lowest order perturbation approximation with respect to 
the effective confining potential V{x)^ if V{x) would be a 
linear function of position over the extent of the unper- 
turbed wavefunction 0^^ yr [x] . The numerical effort with 
this QHA is much less than that required for the full 
Hartree approximation, since no numerical calculation 
of energy eigenvalues and -functions is necessary. Den- 
sity profiles calculated within this QHA are also shown 
in Fig. n as dotted lines. It is seen that the results are 
very similar to those of the full Hartree calculation, in 
particular also the results for the smearing-out of the in- 
compressible strips. Apparently the smearing effect of 
the QHA is even stronger than that of the full Hartree 
approximation. This is understandable, since the Hartree 
wavefunctions are asymmetrically squeezed in space re- 
gions of a rapid variation of V{x)^ and therefore have 
a smaller spatial extent than the unperturbed Landau 
wavefunctions. 

In Fig. we compare the widths of incompressible 
strips with v{x) — 2 for several approximations. The 
line labeled CSG is the analyti cal half - plane result o f 
Chklovskii et al.^, a^.^ — 16\/aodo/^\/^Jc/[t'o^ — ^'.^cIj 
with do ^ 150 nm a depletion length, Vic (= 2) the filling 
factor of the incompressible strip, and vq = 2Ep/hiiJc the 



effective filling factor. ^^ This result agrees well with the 
corresponding self-consistent result of Ref.lllland for suf- 
ficiently low B values (note the inverted B scale in Fig.El) 
also with our present TFPA result for samples of finite 
width. For small B values, the width decreases propor- 
tional to B, and remains finite throughout the figure. 

The result calculated from our QHA is also included 
in Fig. As for the TFPA, we have determined the 
width of the incompressible strips by a simple extrapo- 
lation, using three points next to a plateau to determine 
a plateau edge. For wide plateaus (large B values), the 
QHA width is only slightly smaller than that calculated 
within the TFA. However, with decreasing B the QHA 
width decreases faster and goes to zero at a relatively 
large magnetic field {hujc/Ep « 0.7), far before the in- 
compressible strips with I'^x) = 4 develop in the center 
of the sample. 

These results require a modification of Fig. [21 Within 
the Hartree-type approximation, the width of the incom- 
pressible strip with filling factor 2 shrinks more rapidly 
with decreasing B and vanishes at Hutc/Ep « 0.7. Be- 
tween this B value and the value hujc/Ep w 0.5 no in- 
compressible strips exist in the sample. At still lower B 
values there is a _B-interval in which only incompressible 
strips with local filling factor 4 can exist. This modifica- 
tion required by the QHA is indicated in Fig. [7| below. 

In view of the following it may be interesting to note 
that the essential effects of the Hartree-type approxima- 
tions can be simulated in a very simple way. If we first 
perform a calculation within the TFPA and then take a 
spatial average, e.g., ^{x) = {2X)^^ J_^dx'iy(x + x'), of 
the filling factor profile ^{x), we will smear-out incom- 
pressible strips of a width less than 2A, while incompress- 
ible strips with a width larger than 2A will survive. With 
A of the order of the magnetic length, we will obtain fill- 
ing factor profiles {^{x) very similar to those obtained in 
the Hartrec approximation. 



III. INCOMPRESSIBLE STRIPS AND 
DISTRIBUTION OF DISSIPATIVE CURRENT 

A. The local conductivity model 

We now describe the effect of an applied current on our 
Hall bar system, following again the approach of Ref. H- 
In the presence of a dissipative current /, the electro- 
chemical potential fJ.*{v) will become position dependent, 
and its gradient E = V/x*/e will drive the current density 
j(r). We assume the linear local relation (Ohm's law) 



j{x) ^a{x)'E{x), a{x) =a{nci{x)), 



(14) 



with a position-dependent conductivity tensor d-{x), 
which has the same form as for a homogeneous sample, 
however with the homogeneous density replaced by the 
local electron density nci{x). Due to the translation in- 
variance in the y direction, which is indicated in Eq. (|14|l . 



and the equation of continuity, the components jx and Ey 
of current density and electric field, respectively, must be 
constant j^ 



Mx) = 0, Ey{x) = E°. 
For the other components one finds 

1 ^n „ , . Pnix) 



jyix) 



pi{x) y- 



Ex{x) 



pi{x) y 



(15) 



(16) 



in terms of the longitudinal component pi 



and the Hall component pH 



Pyoi 



Pxx — Pyy 
of the re- 



sistivity tensor p = a . For a given applied current 
/ — J_^dxjy{x) one obtains for the constant electric 
field component along the Hall bar 



E'y=I 



dx 



1 



Pi{x) 



and for the Hall voltage across the sample 

'•' Ph[x) 



/a na 

dxEx{x)=E° / dx 
-d J~d 



Pl{x) 



(17) 



(18) 



With the usual normalization of the resistances to a 
square-shaped conductor, this yields for the Hall and the 
longitudinal resistance 



Rh 



Vh_ 
I 



Ri 



2dEl 



(19) 



Here we see the essence of the local model. Any rea- 
sonable model for the conductivity of a high-mobility 2D 
ES at zero temperature will give simple results for the 
conductivity components at even-integer filling factors 
(where no elastic scattering is possible): 



ai{v=2k) = pi{v=2k) = , 

GH{v=2k) 



PH{v=2k) h 



■2k. 



(20) 
(21) 



Thus, if an incompressible strip of finite width exists in 
the sample, the integral in Eq. (|17|l diverges and Ey and, 
therefore, the longitudinal resistance Ri is zero. At low 
temperatures, ksT <^ hujc, Pi{v=2k) and, therefore, Ri 
will be exponentially small and relevant contributions to 
the integral come only from the incompressible regions. 
The integral in Eq. (|18|l has the same type of singu- 
larity. If only incompressible strips with the same value 
v{x) — 2k of the local filling factor exist, this singular 
integral is just the 2A:-fold of the integral in Eq. (|17|) . so 
that we get the quantized result Rh ~ h/{2ke'^). At zero 
temperature one can evaluate the singular integrals by 
first introducing a cutoff, e.g. by replacing pi{x) with 
pI{x) — max[e,p/(a:)], then calculating the integrals, and 
finally removing the cutoff (e -^ 0). This yields exact 
quantization of the Hall resistance, and the correspond- 
ing calculation at very low temperatures yields exponen- 
tially small corrections. If incompressible strips of finite 



widths with different values of i'{x) exist, e.g. due to a 
manipulation of the background potential, other values 
for the Hall resistance may be possible. But, as we have 
learned from the Hartree-type approximations in the pre- 
vious section, such a situation will not occur in our sim- 
ple translation-invariant Hall bar geometry. From these 
arguments we expect in the resistance-versus-B curves 
plateau regions of finite widths, where the resistances 
have the well known quantized values. 

These considerations are quite general and do not de- 
pend on details of the conductivity model. On the other 
hand, if we want to calculate the resistances between the 
plateau regions, we need to specify a conductivity model, 
and the results will depend on details of this model. We 
will present such detailed results in Sec. II VI below. 



B. Limitations of the local model 

In Sec. ^we have shown that, within a Hartree-type 
approximation, incompressible strips of a width smaller 
than the extent of typical wavefunctions cannot exist. As 
a consequence, incompressible strips with a given filling 
factor 2k do exist only in a finite interval of magnetic 
fields. For lower B values, the i/(x) profile crosses the 
value 2k — v{x2k) with finite slope at some point x ~ X2k- 
At zero temperature, the integrals in Eqs. (|17|) and 118|) 
become singular since pi(x2k) ~ 0. Whether the singu- 
larity is integrable or not depends on the filling-factor- 
dependence of the longitudinal conductivity, (7i{i>). For 
the SCBA model, to be considered below, it is integrable, 
for the Gaussian model considered by GG it is not. But 
should we worry about such sophisticated questions de- 
pending on details of the conductivity model? We think 
we should not, for the following reasons. 

All quantities that are related by Eq. H14|) , the current 
density, the conductivity and the gradient of the elec- 
trochemical potential, represent local values of physical 
variables, which are defined by macroscopic statistical 
arguments. In principle, they have to be calculated as 
average values over sufficiently small subsystems, which 
nevertheless should contain many electrons. We can not 
expect that the local relation (|14|l still holds on a length 
scale of the order of the mean distance between the con- 
stituents of our 2D ES, or, equivalently, of the order of the 
Fermi wavelength Xp- On such a length scale one should 
consider a non-local version of Ohm's law instead. This 
would, however, make things much more complicated, 
and we will not enter such problems. 

In order to simulate qualitatively the expected effects 
of a non-local treatment, we start as before with a local 
model for the conductivity tensor, take the spatial av- 
erage over a length scale of the order of Xp, e.g. with 
X = Xp/2 as 



1 



<?(a^) = ^/ dC&ix + O, 



(22) 



with the averaged conductivity tensor (|22|l . As a conse- 
quence, the resistivity components occurring in Sec. lIH A1 
have to be calculated by tensor inversion of ct(x). 

This simple simulation of non-local effects has several 
appealing aspects. First, if ai{x) vanishes at an isolated 
position X — X2k, the averaged ai{x) > will be pos- 
itive in the neighborhood of X2k, and the integrals in 
Eqs. 117|) and (|18f) will not be singular. Intervals of van- 
ishing ai (x) will exist only, if we start before averaging 
with sufficiently wide incompressible strips (wider than 
2A). Finally, for high- mobility systems, the Hall conduc- 
tivity is given to a very good approximation by the free- 
electron value anix) = {e^ /h)v(x). Thus, the averaged 
Hall conductivity anix) will be given by the averaged 
filling factor profile D{x). As mentioned at the end of 
Sec.^ this averaged profile will agree qualitatively with 
the Hartree profile, if we start with the TFPA profile i'{x) 
and average over X ^ I (l the magnetic length) . Since for 
the large magnetic fields of our interest I '^ Xp, there is 
actually no need to perform the time consuming Hartree 
calculation, if we finally want to calculate the averaged 
conductivity tensor 1)22(1 . 

To summarize, our approximation scheme that simu- 
lates both, the effect of finite width of the wavefunctions 
in the thermal equilibrium calculation, and non-local ef- 
fects on the transport, is as follows. We calculate the 
density profile ^{x) within the self-consistent TFPA, and 
with this the conductivity tensor (t(i/(x)). Then we per- 
form the averaging of Eq. (|22|l and follow the calculations 
described in Sec. lIII Al In contrast to Ref. i^we restrict 
our calculations here to the linear response regime and 
do not investigate the feedback of the applied current 
on the electron density and the electrostatic potential, 
that is mediated in principle by the position-dependent 
electrochemical potential /i*(r) in the presence of a dis- 
sipative current. 



C. Self-consistent Born approximation 

In principle we could use the conductivity models of 
Ref. 5 in order to calculated explicit examples. We prefer, 
however, to take the transport coefficients from the self- 
consistent Born approximation (SCBA),^^'^^'^^ which al- 
lows for consistent models of longitudinal and Hall con- 
ductivities, and for the consideration of anisotropic scat- 
tering by randomly distributed finite-range impurity po- 
tentials. We assume that the relevant scatterers are 
charged donors distributed randomly in a plane paral- 
lel to that of the 2D ES, with an areal density n/, and 
we approximate the impurity potentials by Gaussian po- 
tentials 



v{t) 



Vq 
7ri?2 



exp 



r 

i?2 



(23) 



and use still the local version H14|) of Ohm's law, but now 



with a range R of the order of the the spacing between 
2D ES and doping layer. 



An important aspect of the SCBA is that, similar 
to the "lowest order cumulant approximation" used in 
Ref. 0, it allows us to treat the transport coefBcients and 
the collision broadening of the LLs in a consistent man- 
ner. The relevant SCBA results for the transport coef- 
ficients and the collision broadening of a homogeneous 
2D ES are summarized in the appendix. Consistency 
with the transport coefficients requires that we replace in 
the self-consistent TFPA calculations the (5-functions of 
the Landau DOS (|lll) by the semi-elliptic spectral func- 
tions (jA.Sp . In addition to the range R, the impurity 
strength Vq and concentration n/ determine these quan- 
tities via the relaxation time tq defined by the energy 
H/tq = mVQm/h^. In strong magnetic fields, this en- 
ergy combines with the cyclotron energy to 



r^ = 4n/FoV(2^/') = (2/7r) hwch/ro, 



(24) 



where F is the width of the LLs in the limit of zero-range 
scattering potentials {R ^ 0). We find it convenient to 
characterize the impurity strength by the dimensionless 
ratio 7 — T/huJc at the strong magnetic field _B = lOT 
and define (for GaAs), therefore, the strength parameter 



7j = [{2njV^m/Trh^)/{17.3meY)] 



1/2 



(25) 
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FIG. 4: Filling factor u{x) ^ h aH{x)/e^ versus position in the 
left half of a symmetric high-mobility {R = 20 nm, 71 — 0.1) 
Hall bar of width d = 1.5 fim, for three values of the mag- 
netic field, hiOc/Ep = 0.6, 0.8, and 0.95, and without (A = 0, 
solid black lines) and with (A — 20 nm, long-dashed green 
lines) averaging according to Eq. 11221 . The insets show the 
plateau regions (incompressible strips) enlarged and include 
in addition two results for larger collision broadening but no 
averaging (A = 0). Other specifications in the text. 

Figure. 0] shows the effect of collision broadening on 
the density profile at strong magnetic fields. The sam- 
ple parameters are d = 1.5 /xm, uq — A ■ lO^^cm^^ and 
b/d = 0.952, resulting in fiei = 3.37 • lO^cm-^ and 



t — kBT/Ep = 0.01, three values of the magnetic field 
{huc/Ep — 0.6, 0.8, and 0.95, corresponding to central 
filling factors i^(0) = 3.33, 2.5, and 2.1, respectively), and 
for three sets of the impurity parameters R and 7/. It is 
seen from Fig. 0] and Table Q] that, for sufficiently small 
collision broadening (small 7/ and, eventually, large R), 
incompressible strips still may exist, that their width de- 
creases, however, with increasing broadening of the LLs. 
Table Q] shows, for several sets of impurity parameters, 
the relative widths 7„ of the lowest LLs and the zero 
field mobilities. Data for the second set (i? = 10 nm, 
7/ — 0.1) are not included in Fig. 01 since they could 
not be distinguished from the traces for the first, high- 
mobility set. From the insets of Fig. 0] it is evident that 
incompressible strips can only survive, if the gap between 
the broadened LLs remains broad enough. A large col- 
lision broadening (low-mobility set No. 4) has a similar 
effect as a spatial averaging (long-dashed lines in Fig. QJ 
and may smear out the incompressible strips. 
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0.34 
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0.5 


0.56 


0.40 


0.35 


3.00 



Ep 



12.02 meV. 



E", 



13.51 meV. Data are shown for 



TABLE I: Relative width 7„ = Tn/hwc of the Landau levels 
n = 0, 1, 2 at hujc/Ep — 0.6, and mobility /is=o at _B = 0, 
T = in units of m'^/Vs, for four sets of model parameters 
i?, 77. 



IV. RESULTS AND DISCUSSION 

A. Effect of spatial averaging 

The effect of spatial averaging, introduced to simulate 
non-local effects on the scale of the Fermi wavelength, is 
illustrated in Fig. \E\ It shows, for a magnetic field value 
corresponding to a central filling factor i/(0) = 4.18, the 
filling factor profile calculated with the SCBA broadened 
DOS, together with the conductivity profiles. Fig. [SJa). 
Here we have assumed a short-range potential (leading 
to the rather low mobility fJ,B=o = 6.4m^/Vs), in or- 
der to obtain a noticeable deviation of the (Jh{x) trace 
from that for the filling factor v{x). Clear incompress- 
ible strips with the quantized values for cjh{x) and van- 
ishing ai{x) are visible where v{x) assumes the integer 
values 4 and 2. The effect of spatial averaging on the 
conductivities is demonstrated in (b) and (c). The wide 
{^ 90 nm) plateau defined by i'{x) = 4 shrinks due to 
the averaging (to ^ 50 nm) but clearly survives, as is 
shown in Fig. |SJb) for cri(x), and hold similarly for the 
plateau of anix). On the other hand, the plateau behav- 
ior near the narrow (~ 25 nm) strip defined by i'(x) = 2 
is completely smeared out, and cri(x) does no longer van- 
ish in this region, Fig.[S{c). This has, of course, drastic 
consequences for the current distribution, which is dom- 
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FIG. 5: (a) Filling factor and conductivity profiles for the 
left half of a symmetric sample with d — 1.5 /j-m, no — 
4 • 10"cm~^, 6 = 0.9d, calculated within the SCBA, for R = 
0.1 nm and 7/ = 0.1, hujc/E°p = 0.48, and ksT/E^ = 0.01. 
(b) and (c) repeat the data of (a) (solid lines) near the in- 
compressible strips with filling factors ^{x) — 4 and ^{x) = 2, 
respectively. The dash-dotted lines demonstrate the effect of 
spatial averaging, according to Eq. 1)22^ , with A = 30 nm. 



inated by the strips with vanishing ai{x), i.e., vanishing 
pi (x) , see Eq. (fT?)| . Without averaging a finite fraction of 
the total current would flow through the incompressible 
strips with i'{x) = 2 (on both sides of the sample). With 
the averaged conductivity tensor, the total current must 
flow through the incompressible strips with iy{x) — 4 (as- 
suming that there ai (x) = holds exactly) . 

The mechanism illustrated in Fig. is, of course, also 
effective at other values of the magnetic field. At larger B 
with i^(0) > 2, broad incompressible strips with 1^(2;) = 2 
will exist near the center of the Hall bar, and the spatially 
averaged conductivities will show clear plateau behavior. 
With decreasing B, the incompressible strips move from 
the center towards the sample edges and shrink. If the 
strip width becomes of the order of 2A or smaller, the 
averaging according to Eq. l|^ will destroy the plateau 
behavior of the conductivities and ai{x) will no longer 
vanish near the strips. Then the current density still may 
have a (finite) maximum near the strips with ^{x) — 2, 
but a finite amount of current will spread over the bulk 
of the sample, and the global resistances will no longer 
have the quantized values Rh = h/2e'^ and Ri = 0. 

Figure shows typical results for the dependences of 
the global resistances on magnetic field and averaging 
length, as calculated from Eq. H19|l . At high magnetic 
fields, hujc > Ep, we have everywhere in the sample 
i'{x) < 2, and no incompressible strips exist. The fill- 
ing factor and consequently the conductivities and the 
current density vary slowly across the sample. Thus, the 
spatial averaging has little effect and, within the accu- 



FIG. 6: Calculated Hall and longitudinal resistances versus 
scaled magnetic field hcuc/Ep, for different values of the av- 
eraging length A. The sample parameters are d — 1.5 /xm, 
no = 4 ■ 10"cm"2, b = 0.952d, _R = lOnm and 77 = 0.1, and 
kBT/E% = 0.01. 



racy of the figure, the results with and without averag- 
ing agree. For slightly lower magnetic fields, hwc ^ Ep, 
broad incompressible strips exist near the center of the 
sample and, for all the considered averaging lengths A, 
strips of finite width with ai{x) = and <th{x) = 2e^//i 
survive. As a consequence, the resistances are quantized, 
independently on the width of these strips. 

For still lower B values the situation becomes more 
complicated. Within the TFPA, incompressible strips 
exist for all these B values. Without spatial averag- 
ing, ai [x) vanishes on these strips and, as a consequence, 
Ri = 0, and, for hwc/E^p > 0.5, Rh = /i/2e^, as shown 
by the traces for A = 0. For hujc/Ep < 0.5 there are two 
types of incompressible strips, with ajiix) = 2e^/h or 
Ae^/h, and, without averaging, both contribute accord- 
ing to their widths to Rh, while still Ri = 0. The fluctu- 
ations in the Rh curve for A = result from our TFPA 
calculation on a finite mesh (A'' = 500), which yields dis- 
continuous changes of the widths of the incompressible 
strips with changing B. This unsatisfactory picture, ob- 
tained for A = 0, results from the model calculation of 
Ref. '5. 

The introduction of the spatial averaging improves this 
situation dramatically and leads to qualitatively correct 
results. With decreasing B the incompressible strips 
with ^{x) = 2 become narrower. As the width becomes 
smaller than 2A, ai{x) no longer vanishes, the integrals 
in Eqs. H17|) and (|18|l and thus Ri become finite. This 
happens at higher B values if A is larger, and the resis- 
tances near the low-magnetic-field edge of the quantum 
Hall plateau depend strongly on A. While ai{x) may 
have a sharp minimum near the strips with ^{x) = 2 if 
the width of these strips is only slightly smaller that 2A, 
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FIG. 7: Calculated Hall resistance (light solid line) and 
(scaled) longitudinal resistance (black solid line) versus mag- 
netic field, measured in units of huc/Ep. The underlying gray 
scale plot shows the averaged filling factor profile, as in Fig.|21 
The crescent-like areas indicate the regions of incompressible 
strips with local filling factors 2 (right) and 4 (left). 



this mininium, and also the corresponding maximum of 
the current density, will smoothen out as the width of the 
strips becomes much smaller than 2A. Then the total re- 
sistances will become nearly independent of A, as is seen 
in Fig. El for 0.5 < hujjE^p < 0.6. For hujjE'^p < 0.5, 
ai (x) vanishes only within the incompressible strips with 
^{x) = 4, but not in strips with ^{x) — 2. As a conse- 
quence, we obtain again the exactly quantized results for 
Rh and Ri. 

To visualize the intimate connection between the ex- 
istence of incompressible strips of finite width [now with 
constant anix) ~ e'^i?{x)/h], we show in Fig. [7| a gray 
scale plot of the spatially averaged filling factor profile 
for a relevant interval of magnetic fields, together with 
the resulting resistances. 



B. Effect of temperature and collision broadening 

The spatial averaging procedure is essential to obtain 
quantum Hall (QH) plateaus of finite width for the Ri{B) 
and the Rh{B) curves and to obtain the correct quan- 
tized values for Rh{B) on the plateaus corresponding to 
filling factors larger than two. The width of the calcu- 
lated QH plateaus does, however, not only depend on the 
averaging length A, but also on the temperature and on 
the broadening of the LLs due to the impurity scattering, 
since both affect the width of the incompressible strips. 
The effect of collision broadening on the width of incom- 
pressible strips has been indicated in Fig. 0] The tem- 
perature effect has been investigated in Ref. [l^, where 
it was shown that, in the absence of collision broaden- 
ing, i.e., on the basis of the bare Landau DOS (|ll|l . the 
incompressible strips have a finite width at zero tem- 
perature. At finite, increasing temperatures, the width 



shrinks (while the value of the filling factor remains ex- 
actly constant within the remaining strip) until at a suffi- 
ciently high temperature {ksT < huJc/25 in Ref. lUi ) the 
width collapses to zero. A similar result is expected in 
the presence of collision broadening. In particular we ex- 
pect that, within the self-consistent TFPA based on the 
SCBA DOS, the existence of an energy gap between two 
adjacent broadened LLs will always lead to an incom- 
pressible strip, provided the temperature is low enough. 
The necessary temperature will be the lower, the nar- 
rowtygt he gap is. Hut we^ -yyill r|ot discuss these .q'lCstin ns 



in further details 
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FIG. 8: Hall and longitudinal resistances versus magnetic 
field, calculated for different temperatures, t = kBT/E%. 
Sample parameters: d = LS/im, no — A ■ 10^ cm~ , 6 = 
0.952d, i? = lOnm and 77 = 0.1, A = 30 nm. 

The temperature effect on the calculated resistance 
curves is shown in Fig. |S1 As expected, the width of 
the QH plateaus increases monotonically with decreas- 
ing temperature, but apparently it has a finite limit for 
T ^- 0. We have also included the high-temperature re- 
sult for kBT/E°p = 0.3 (since E°p = 13.5 meV, this means 
T « 47 K). At the low-B side of the figure ksT - hiUc, 
and the derivative of the Fermi function overlaps about 
two LLs. Since we consider only the lowest Landau lev- 
els (n = 0, 1, 2), our calculation is not correct in this 
limit. Nevertheless it is interesting to compare this re- 
sult with the Drude result, which should be valid if the 
Shubnikov-de Haas oscillations are smeared out at higher 
the temperature. 

In the Drude approximation we have Ph{x) = 
^cTtrPi{x), with pi{x) = l/cro(x), where 0-0(2;) = 
e^ Ttrnc\{x) / m is independent of the magnetic field. In- 
serting this into Eq. (|17|l . we obtain 



/ 



e^Ttr e2 2Ef 

£/,'' m n n/Ttr 



(26) 



In Eq. (|18|l we take the integrand to be tOcTtr, but only 
for l^l < b, where nei(a;) is not exponentially small, and 
obtain Vh/E^ = 2buJcTtr- With Eq. (^3 we obtain the 
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Drude result 









(27) 



The energies Ep and E'p are calculated numerically from 
the density profile at B = 0, T = 0, as is described above, 
and Ttr is calculated as described in the appendix, with 
kp = ^/27rnoi. The results are plotted as dashed straight 
lines in Fig. |S1 

Finally, Fig. |2| shows the effect of the Landau level 
broadening on the QH plateaus at a fixed, relatively low 
temperature. The corresponding widths of the lowest 
LLs, in units of the cyclotron energy fiwc, are given in 
Table ^ for the lowest and the largest B value shown 
in the figure. For the largest damping, the LLs start to 
overlap for hwc/E^ < 0.4. so that the v = A QH plateau 
is not well developed. We note that the SCBA results 
summarized in the appendix are only valid at higher B 
values, where the LLs do not overlap. 

The two high-mobility situation considered in Fig. 
differ only by the range R of the Gaussian impurity po- 
tentials. The larger range leads to slightly smaller level 
broadening, but to much lower longitudinal resistance 
(i.e., to much higher mobility at i? = 0). 
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FIG. 9: Hall and longitudinal resistances versus magnetic 
field, calculated for different values of the collision broadening. 
Sample parameters: d = 1.5 fim, no — i ■ lO^^cm"^, b — 
0.952d, A = 30 nm, kBT/E"p = 0.01 (i.e., T = 1.57K). 



C. Hall potential profile 

The motivation of Ref. |a and our present work came 
from the experimental investigation^ of the electrostatic 
potential distribution across a Hall bar under QH con- 
ditions, caused by an applied current. Ahlswede and 
coworkersiii observed three types of potential distribu- 
tion, depending on the filling factor regime. Type I was 
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TABLE II: Relative widths of the lowest Landau levels, j^ = 
r„/riujc for hu}c/E% = 0.35, and 7^ = r„/hujc for hwc/E% = 
1.10, for the impurity parameters used in Fig.|^ First column: 
range in nm, last column: zero field mobility in ra? /Vs. 



a more or less linear variation across the sample and is 
observed if the filling factor in the center is smaller and 
relatively close to (but not too close to) an integer n, i.e., 
n > iy{0) > n — 1/2. If the center filling is slightly larger 
than an integer, n < z/(0) < n-|- 1/2, type HI is observed, 
characterized by a constant potential in the central re- 
gion and a rapid variation across (narrow) strips, which 
move with decreasing B towards the sample edges and 
have been interpreted as incompressible strips^ Finally, 
type II shows a rapid, non-linear variation of the po- 
tential in the center region and is observed if the center 
filling factor is very close to an integer. 

In Ref. la it was shown that, in an local equilibrium 
picture, the changes of the electrostatic potential, caused 
by an applied current, follows closely the current-induced 
variation of the electro-chemical potential n* , so that the 
resulting density changes are small. In the present work 
we do not consider the feed-back of the spatial variation 
of fi* on electrostatic potential and density profiles(linear 
response). But we expect from the results of Ref.iSl, that, 
in the linear response regime, fi* should show the same 
position dependence as the electrostatic potential would 
do, if the feed-back were calculated. 

To calculate the Hall profile across the sample, we in- 
tegrate Ex{x), Eq. H16|l . from the center x = to the 
actual X value. Typical results as functions of position x 
and magnetic field B are shown in Fig. 1101 The normal- 
ization is chosen so that Vh{B,x = d)/I = —Vh{B,x = 
—d)/I = Rh/'^- One sees clearly that the plateaus 
of the quantized Hall effect (0.8 < hLJc/Ep < 1 and 
0.45 < hujc/Ep < 0.5) coincide with potential variation 
of type III, caused by current density confinement to the 
incompressible strips. Moving from a plateau region to 
smaller B values, the incompressible strips shrink and fi- 
nally vanish, and the current density spreads more and 
more out into the bulk. This leads to the type I behavior 
(0.52 < hLJc/E"p < 0.7 and hwc/E^ < 0.4). Immediately 
above the integer values of the center filling factor (in our 
approximation assuming spin degeneracy near ^{O) — 2 
and 4), we find the rapid variation of the type II. This 
is in very nice agreement with the experiment. With- 
out our spatial averaging of the conductivity tensor, we 
would have missed the type I regions for hujc/Ep < 1, as 
has been observed in Ref. 15|- 
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FIG. 10: Hall potential profile Vnix) = J^ dx'Ex{x') across 
the sample, for varying B/Bq = hujc/ E% and constant applied 
current /. Normalization: \VH[d) — Vh(— d)]/J = Rh\ sample 
parameters: d = 1.5 //m, no = 4 ■ lO^^cm" , h = 0.952d, 
A = 30 nm, ksT/E^ = 0.01. 



D. Summary 

The virtue of our approach is, that it allows us to 
calculate resistance traces with exactly quantized quan- 
tum Hall plateaus of finite width, and with reasonable 
values of the resistances between these plateaus. While 
these intermediate resistance values depend on the details 
of our conductivity model, the quantized plateau values 
do not. The reason for this high accuracy and model- 
independence of the plateau values is the fact, that the 
latter are determined by the integrals in Eqs. H17() and 
(|18|l becoming singular across incompressible strips. 

To obtain realistic widths of the QH plateaus, we had 
to consider a mechanism that prohibits singular current 
flow along very narrow incompressible strips. We have ar- 
gued that small-scale non-local transport effects act into 
this direction, and that consideration of the finite extent 
of wavefunctions will prohibit arbitrarily narrow ISs at 
low magnetic fields, in contrast to the prediction of the 
Thomas-Fermi approximation. We were able to simu- 
late such non-local effects by a simple spatial averaging 
procedure, with reasonable results for Hall and longitu- 
dinal resistance as functions of the magnetic field. Also 
the resulting potential profile, and therefore the current 
distribution across the sample, is in nice agreement with 
recent investigations.^ We consider this as a strong sup- 
port for the relevance of our approach, notably because 
earlier approaches, which neglected dissipation, cannot 
explain the experiments, as has been discussed in Ref. 'S". 

Note that, for QH plateaus corresponding to filling fac- 



tors v > i, our results are qualitatively different form the 
conventional edge channel picture. The latter explains, 
for instance, the quantized conductance value G = 4 e^/h 
as the sum of the contributions of two spin-degenerate, 
quasi-one-dimensional current channels near each of the 
opposite sample edges, thus tracing back the quantized 
Hall effect to the phenomenon of ID conductance quanti- 
zation (in a situation where no backscattering occurs)^ 
That is, the edge states, created by the LLs with quan- 
tum numbers n = and n = 1, contribute both to the 
current in the plateau regime of the QHE. Our results, on 
the contrary, indicate that the total current flows along 
the incompressible strip with local filling factor i/^x) = 4 
(where both LLs n — and n — 1 are occupied) , whereas 
near local filling factor i^(x) = 2 no incompressible strip 
and no contribution to the current exists. 

Comparing our resistance curves with experiments, we 
notice that the high-field edge of a calculated plateau oc- 
curs at a magnetic field, at which an incompressible strip 
(with an even integer value of the effective filling factor 
vq = 2Ep/hiOc) first occurs in the center of the sample. In 
experiments these i^q values usually are found somewhere 
near the centers of the plateaus. We have good arguments 
that this discrepancy is due to our neglect of long-range 
potential fluctuations due to the randomly distributed 
ionized donors. We have simulated the "short-range" 
part of the Coulomb potentials of the remote donors by 
Gaussian potentials, but we have neglected their overlap- 
ping long-range parts, which lead to long-range potential 
fluctuations. We have evaluated the short-range disor- 
der within the SCBA to calculate conductivities and LL 
broadening. We have seen that with increasing disor- 
der scattering the level broadening increases and, as a 
result, the widths of the QH plateaus shrinks. On the 
other hand, one knows from technical applications of the 
QH effect, that rather impure samples have usually es- 
pecially wide and stable QH plateaus. This points to the 
role of long-range potential fluctuations, which become 
more important with increasing impurity concentration. 

As a rough simulation of such long-range fluctuations, 
we have added oscillatory terms to the confinement po- 
tential and than repeated our calculations. We indeed 
find that such modulations can widen and stabilize the 
QH plateaus, and eventually even shift them to higher 
magnetic fields, depending on the amplitude, the range, 
and possible other details of the perturbation. This be- 
comes understandable, if one considers the effect of such 
fluctuations on the existence and the position of incom- 
pressible strips. However, we do not want to discuss such 
considerations further, since we believe that our quasi- 
one-dimensional model is not appropriate for a reliable 
discussion of statistical fluctuations of a 2D donor dis- 
tribution. Effects of long-range disorder in an uncon- 
fined 2D ES on the longitudinal resistance between QH 
plateaus have already been discussed a decade agoii^ 
This discussion seems however not applicable to the 
rather narrow samples of our present interest, since, first, 
the confinement affects the self-consistently calculated 
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potential and thus the density distribution,— and, sec- 
ond, the early assumptions about current-carrying and 
isolating regions are not compatible with our results and 
the experimental findings^ 
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APPENDIX: SCBA CONDUCTIVITIES 



The low-temperature, high-field magnetotransport, de- 
termined by elastic scattering of the 2D electrons by 
randomly distributed impurities with scattering poten- 
tials of arbitrary range, has been studied by Ando 
and coworkers i^^'^^i^^ The results for the case of non- 
overlapping LLs can be summarized as 




FIG. 11: SCBA results for longitudinal [ai) and Hall (an) 
conductivity, in units of e^/ft, versus filling factor, at fixed 
magnetic field for different values of impurity range (a = 
R/l) and strength (7 — T/huc), and of temperature (t = 
ksT/huJc)- For the two lower panels the correction —Agh to 
a% = e^f/h is negligible, and an is not shown. 
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Assuming a single type of impurities with the Gaus- 
sian potential H23|) . these parameters can be expressed 
in terms of the integrals 



(r(/))^ 



dx g^-'') (x) cxp(- [1 + a^]x) , (A.8) 



where F^ = AniVQ/{2TTp) and a = R/l, and the weight 
functions 






2.^ ^^"^(-)' 



are determined by the associated Laguerre polynomials 
L™{x). With these notations one obtains 

f2=(f(«))2, (F--)2 = (FW)^ (A.9) 

(Fr)^ = (FWr + (F(-))^ (A.IO) 



In the limit of short-range scattering potentials, a ^ 0, 
one has F^/F^ = 1, (F^^/F)^ = n + l/2 and (F^^^/F)" = 
ri+1/2. With increasing a these parameters decrease and 
remain for a > 1 about an order of magnitude smaller 
than their a = values. ^"^'^^ 

Some typical SCBA results are shown in Fig. ^2 In 
general, cr; and the correction Act// to the free electron 
Hall conductivity a'j^ = e^v/h decrease with increasing 
range of the scattering potentials. 

At zero temperature, Act/t is proportional to T/huJc 
with a factor depending only on the range parame- 
ter a = R/l. The longitudinal conductivity (Ji{v), on 
the other hand, depends only on a and not on F^ = 
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(2/Tr)HLOcfi/To, i.e. depends not on the impurity concen- 
tration 71/ and strength Vq entering the B — relaxation 
rate h/ro — niV^^m/h^. This is very different from the 
B — Q conductivity ctq — e^n^xTtrl'm obtained for the 
impurity model H23I) . which depends on rid = k^j'l'K and 
via Ttr on both the potential strength and range. For 
i3 = collision broadening effects can be neglected, and 
one obtains for elastic scattering at the Fermi egde 

% njm f^ dip 



Ttr 



m r dip r -.2 , . 

— y_^^Nq=fe^[2(l-co«^)]i/2(l-COS(p) 



[e-^I,{x)-h{x))]^__ 



(A.ll) 



where the last equality holds for our impurity potential 



(|23|l . with Fourier transform Vq = Vq exp(— i?^(7^/4), and 
Ii, (x) is a modified Bessel function. This leads for Rkp ^ 
1 to To/rtr ~ [\/8n (Rkp)^]^^ ■ With increasing temper- 
ature, the peak values of ai{i') decrease and the minima 
at even integer v values are no longer exponentially small 
for ksT/hwc > 0.1. This behavior of the SCBA results, 
shown in the middle panel of Fig.^1 is similar to that of 
the Gaussian model treated in Ref . ^. At finite tempera- 
ture, the longitudinal conductivities cri{h') increase with 
7, i.e., with increasing scattering strength (bottom panel 
of Fig. Ill|l . This is as expected from the Drude picture 

for LOcTtr > 1: (Ti W ao/{LLlcTtr)'^ CX l/Tfr- 
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